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Abstract 

We give an explicit solution describing internal waves with a still wa- 
ter surface, modelling the dead water phenomenon, on the basis of the 
Gerstner wave solution to the Euler equations. 

1 Introduction 

The phenomenon of dead water was first investigated by Vagn Walfrid Ekman 
at the initiative of the Norwegian oceanographer and explorer Fridtjof Nansen, 
who encountered it on the voyage of the Fram through the Arctic Ocean, and 
subsequently proposed its study to Ekman's teacher Vilhelm Bjerknes. Descrip- 
tions of dead water as encountered by early seafarers go back to antiquity, an 
engaging ac count of whic h is given by Ekman in The Norwegian North Polar 
Expedition ( Nansen . 1906t h 



Dead water is so-called for its ability to slacken the speed of a vessel quite 
suddenly, and with no apparent cause on the water surface. This phenomenon 
is prevalent when a layer of less dense water overlays a layer of denser water, 
due either to differences in salinity or temperature. Such instances of fresh or 
brackish water resting upon heavier sea- water occur at the mouths of rivers, 
particularly in the Norwegian fjords, where relatively little mixing of the waters 
takes place. Internal waves generated at the interface between the two fluids 
may impede the progress and steerage of ships. 

Causes and countermeasures for dead water reported by early mariners were 
manifold and at times quite amusing. It was believed that the dead water could 
be loosened from a ship by beating the water with oars or pouring petroleum 
into the water ahead of the ship. The belief that it was caused by a small remora 
fish sucking fast to the ship's body was also widespread, an explanation tracing 
back to Pliny the Naturalist, and there abound stories of submarine magnetic 
rocks, or supernatural forces. 

In the tradition of the early 20" 1 century, Ekman's work r elied on a linea rized 
theory (based on the discussion of internal wave s by G . G. Stokes ( 18471 )'). as 



did Lamb's article on the subject a decade later ( Lamb . 1916t . More recently, 



however, both theoretical and experimental investigations have highlighted the 
importance of nonlinear effects in connection with the dead water phenomenon. 
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Mercier, Vasseur and Dauxois' experimental study ([Mercier et al I l201lh revisits 
Ekman's work and goes on to study interfacial waves in three-layer as well as 
linearly st ratified fl u ids w ith a pycnocline. Motivated in part by this experimen- 
tal study, iDuchene ( 2011 ) has given a sophisticated analysis of some nonlinear 
models for the dead-water phenomenon with the assumption that the wave- 
lengths of the internal waves are long compared to the layer depths, while also 
incorporating aspects of the ship motion. We refer the reader also to references 
therein for a fuller account of the existing literature. 

We will provide an explicit solution to the nonlinear governing equations 
for water waves describing such internal waves which leave no trace upon the 
surface of the water. While the passage of a ship will necessarily generate waves 
both on the water surface and at th e pycno c line, the surface waves will be 
of considerably diminished amplitude iNanseii ( 1906 ). and these long and low 
waves may we l l esca pe notice both by mariners as well as in the laboratory 
( Mercier et alll201ll p. 194). 



The mechanism of generation of these waves by a vessel is, regrettably, out- 
side of the scope of this study. Hence we must content ourselves with a descrip- 
tion of the wave motion itself, based on the simple insight that the classical 
Gerstner wave can propagate also at an interface between two liquid regions. 
The difference between the problem presented at the two-fluid interface and 
the usual water-wave problem at the free surface manifests in the fact that the 
pressure at the interface is no longer constant, but hydrostatic. The Gerstner 
wave construction allows the pressures at the interface to be matched, which 
introduces a reduced gravity into the problem. 



2 Governing equations 

The physical situation as seen in figure 1 may be captured as follows: Let our 
fluid be contained in the region {(x,y) | x £ R, y < 0}, while the upper half- 
plane y > is assumed to be composed of air with negligible density which docs 
not interact with the fluid below. We will assume that the upper layer of fluid 
has a constant density pi, and the lower fluid a density p 2 , with p\ < p 2l the 
interface between these immiscible layers being denoted by rj(x,t). 




Pi 



Figure 1: Internal waves at the interface between lighter (density pi) and heavier 
(density p 2 ) water. 
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The governing equations for two-dimensional, incompressible, inviscid free 
surface flow are the equation of mass conservation 



(1) 



and the Euler equation 



p p 

u t +uu x + vu y = — v t + uv x + vv y = -— - g, (2) 



supplemented by suitable boundary conditions cf. ( Constantmll2012h . Here the 



horizontal and vertical components of the velocity field are u(x, y, t) respectively 
v(x,y,t), while P(x,y,t) denotes the pressure, and g is the constant gravita- 
tional acceleration. Subscripts indicate a partial derivative. 

In order to decouple the motion of the air from that of the water, we intro- 
duce a dynamic boundary condition at the water surface: 

P = P atm on y = 0, (3) 

which specifics that the pressure at the upper interface is equal to the atmo- 
spheric pressure P a tm- Wc are interested in an interface describing steady waves, 
and may thus assume rj = rj{x + ct). In order to assure that this interface sepa- 
rates the fluids completely, one introduces a further, kinematic, boundary condi- 
tion specifying that the fluid velocity thereon is wholly tangential. Temporarily 
denoting by (u 1 ^ 1 ) and (u 2 ,v 2 ) the velocity fields of the upper respectively 
lower layer, this condition takes the form 

v l = ?/(«' + c) on y = rj, for i = 1, 2. (4) 

As we shall not account for motion far below the water surface, there is no need 
to specify a similar condition for the sea-floor. Instead, we will assume that our 
motion dies out at great depth 

u 2 , v 2 — » as y — > — oo. 

The top layer of water will be assumed to be moving with the wave at speed 
c and the surface of the water taken to be flat, such that (u 1 , v 1 ) = (— c, 0). This 
satisfies the kinematic boundary condition at the interface (U), and, trivially, 
the same condition at the water surface. Hence, in the near-surface layer, ([2]) 
becomes 

P x = and — P y = pig 

which, along with the condition ([3|) that the pressure be equal to the atmospheric 
pressure at the flat water surface, yields 

P = Patm - pigy (5) 
throughout the upper fluid layer. 
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3 Structure of the solution 



In order to describe the internal wave we are interested in, we specify the par- 
ticle paths of our solution in Lagrangian, or material, coordinates (a, 6), with 
a € K and b < bo < 0, where bo will describe the interface between the two 
fluids. Indeed, t he wave motion is simply that of Gerstner's trochoidal solution 
(|Gerstneii [l809l) . though we shall see that the pressure at the interface 77 will 
lead to some modifications. The particle trajectories are given by 

e mb . 

x = a H sinmfa + ct), 

TO 

e mb 

y = b cos m(a + ct). 

TO 

Then the horizontal and vertical velocities are 

u = ce mb cos m(a + ct), 
v = ce mb sin to (a + ct ) , 

and we find the accelerations to be: 

= —c 2 me mb sin m(a + ct), 
Dt y ' 

-j^- = c me cos m{a + ct). 

Because the motion is two-dimensional, we may also specify the velocity field 
via a stream function, 



tp = c ( b — 



„2mb 

2m 

whereupon we may identify the streamlines to be curves of constant b. Setting 
for readability to (a + ct) = 9 and mb = £, we find 

Px = P2c 2 me^ sin 9, 

P y = —p2 (c 2 me^ cos 6 + g) . 

Now, transforming to material coordinates via 



(6) 



P a \ _ /l + e^cosfl e« sin 9 \ (P x 
P b )~\ e«sin6» 1 - e« cos 9) \P y 

yields 

P a = P2 (c 2 m - g)e^ sin 9, 
Pb = P2<?me 2 ^ — P2.9 + P2{g — c 2 m)e^ cos 9, 
which we may integrate to yield a pressure 

P = -p 2 (c 2 m- g)—cos9 + p 2 —e 2 t - p 2 gb + C. (7) 

TO 2 
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At the interface, which we take as the streamline b = bo, the above pressure 
must match that for the still fluid ([5]), given in material coordinates by 

/ e mb \ 
P = Patm - pig ( b - — cos m(a + ct)j . (8) 

This means that we must require 

9 Pi = P2(g - c 2 m), 

which is equivalent to 

2 02 - Pi 
c m = g , 

P2 

the right-hand side of which is simply the reduced gravity, which we denote by 
go, implying that the interfacial waves propagate with a celerity 

c = \/go/m. 

This is clearly consistent with the propagation speed of classical Gerstner waves, 
which is recovered by setting p\ = 0. The equation of mass conservation ([!} 
is equivalent to the determinant of the coordinate transform ([5]) being time 
independent, which is readily verified. 

The parameter bo which defines the interface also determines the form thereof: 
6o = corresponding to a cycloid with sharp crests, while bo < describes a 
smooth, trochoidal profile. The assumption that bo < is made to ensure that 
there are no self-intersecting particle paths. The matching of pressures (|7|) and 
([8j at b — bo further necessitates that bo be the solution to the equation 

^-b)=P atm -C. 

This solution is unique for a suitable choice of constants, as the left-hand side 
is a strictly decreasing function of b whose value in the limit b — > — oo is clearly 
+oo. 




4 Discussion 



We have provided a solution describing waves with a trochoidal profile prop- 
agating at the interface between water of two different densities with a flat 
surface. The Gerstner wave which forms the cornerstone of this study of dead- 
water belongs to t he early history of water-wave theory, first discovered by its 
namesake in 1804 (IGerstnerl . 11809) . and rediscov ered by William J. M. Rankine 



more than a half-century later Plankinel . Il863h . A modern treatment of the 
Ger stner fl o w, pro ving its dynamic feasibility, was given bv lConstantinl (<200lh 
and iHerirv ( 2008a ). The particle motion is along circular trajectories with an- 
gular velocity mc, the radius decreasing exponentially with depth. For positive 
celerity c the waves move in the negative x direction, and there is a negative 
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vorticity associated with the wave motion - this latter indicating that the waves 
cannot be generated from rest by potential forces in an inviscid, homogeneous 
fluid. Moreover, the circular particle paths are a hallmark of vorticity, since for 
irrotational Stokes waves the particle trajectories have quite a diffe rent pattern 
cf. rece nt work by IConstanti rJ (|2006h . lConstantin fc Strauss! (|2010l ). and iHenrvl 
( 200861) . 

It should be noted that the peculiar structure of the Gerstner solution is, 
perhaps, a mixed blessing - while it furnishes us with an explicit solution whose 
profile is more realistic than the sinisoid of the linear theory, at the same time 
the nonlinear boundary conditions and governing equations are very exacting, 
leaving only little leeway to include other effects. At the same time, as the Ger- 
stner wave is the only available explicit solution to the full water-wave problem, 
it has been modified to describe a variety of physic al situations , (see f or example 
recent work by Constantin on geophysical waves IConstantin (2013) and refer- 
en ces therein). Desp ite this special structure, recent work on deep-water waves 
by iMonismith ( 2007 ) et al finds that Gerstner wave theory fits observed mean 
velocities better than classical deep-water Stokes theory, both in the laboratory 
and in the ocean. 

In contrast to the classical Gerstner wave solution, with the inclusion of an 
upper layer of still water the pressure ([7]) is no longer constant along streamlines, 
but exhibits an oscillatory time-dependence. On account of this, our solution is 
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also distinc t from the Gerstner wave in stratified water ([Dubreil-Jacotin 
Stuhlmeier . 201ll) . where lines of constant density and constant pressure coin- 
cide with the streamlines. In fact, upon passing to a reference frame moving 
with speed c, constancy of the density along streamlines, i.e. p a = 0, is ensured 
by the continuity equation - in our setting, the form of the pressure precludes 
the possibility of exploiting this to describe a heterogeneous bottom layer. 
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